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ABSTRACT 

The propeller regime of disk accretion to a rapidly rotating magnetized star is investigated here for the 
first time by axisymmetric 2.5D magnetohydrodynamic simulations. An expanded, closed magnetosphere 
forms in which the magnetic field is predominantly toroidal. A smaller fraction of the star's poloidal 
magnetic flux inflates vertically, forming a magnetically dominated tower. Matter accumulates in the 
equatorial region outside magnetosphere and accretes to the star quasi-periodically through elongated 
funnel streams which cause the magnetic field to reconnect. The star spins-down owing to the interaction 
of the closed magnetosphere with the disk. For the considered conditions, the spin-down torque varies 



with the angular velocity of the star (w*) as 



a. 3 



for fixed mass accretion rate. The propeller stage 



may be important in the evolution of X-ray pulsars, cataclysmic variables and young stars. In particular, 
it may explain the present slow rotation of the classical T Tauri stars. 

Subject headings: accretion, dipole — plasmas — magnetic fields — stars: magnetic fields — X-rays: 
stars 



1. INTRODUCTION 

For a rapidly rotating magnetized star, the centrifugal 
force at the magnetospheric radius may be larger than the 
gravitational force. In this case the magnetosphere acts 
as a "propeller" pushing away incoming matter (David- 
son & Ostriker 1973; Illarionov & Sunyaev 1975; Lipunov 
1992). The propeller regime has been discussed in con- 
nection with different astrophysical situations (e.g., Stella, 
White, & Rosner 1986; Treves, Colpi & Lipunov 1993; Cui 
1997; Alpar 2001; Eksi & Alpar 2003; Mori & Ruderman 
2003). 

A detailed understanding of the propeller regime is im- 
portant for the analysis of a number of astrophysical sys- 
tems, including Classical T Tauri stars (CTTS), cata- 
clysmic variables, X-ray pulsars, isolated neutron stars, 
and possibly anomalous X-ray pulsars. 

For quasi-spherical accretion, the propeller regime was 
investigated analytically by Davies, Fabian, and Pringle 
(1979), Davies and Pringle (1981), and Ikhsanov (2002), 
and numerically by Romanova et al. (2003a) for the case 
of an aligned magnetic rotator. The numerical simulations 
have shown that the propelling star expels significant mag- 
netic flux in the radial direction, thus forming the low- 
density, azimuthally wrapped equatorial magnetic wind. 
At the same time, some matter was observed to accrete 
along the magnetic poles of the star (as observed in the 
non-rotating case, Toropina et al. 2003). These simula- 
tions revealed a complicated flow geometry not considered 
in the theoretical models. 

Disk accretion to a star in the propeller regime was in- 
vestigated theoretically by Li and Wickramasinghe (1997), 
Lovelace, Romanova, and Bisnovatyi-Kogan (1999). Wang 



and Robertson (1985) analyzed numerically the disk-star 
interaction in the propeller regime, but only in the equa- 
torial plane. Romanova et al. (2002; hereafter RUKL) 
did axisymmetric simulations of disk accretion to the fast 
rotating star, however only in the regime of very soft pro- 
peller. In this paper we discuss for the first time ax- 
isymmetric MHD simulations of disk accretion to fast pro- 
pellers. 

2. THE MODEL 

Axisymmetric MHD simulations have been done for the 
interaction of an accretion disk with a rapidly rotating 
magnetized star in the propeller regime. The arrangement 
of the model is similar to that used in RUKL for analy- 
sis of disk accretion to slowly rotating magnetized stars. 
That is, (1) a spherical grid (R,9,<fr) was used which as- 
sures the high spatial resolution near the star; the grid 
Nr x Ng — 131 x 51 was used in most of simulations; (2) a 
Godunov-type numerical method was used; the code was 
improved compared to RUKL so that now we were able 
to deal with extremely fast rotating magnetized stars (see 
details in Ustyugova et al. 2004); (3) "quiescent" initial 
conditions were used so that we were able to investigate 
the slow viscous accretion of disk matter to the star. 

We use the same scaling parameters as in RUKL. The 
length-scale Rq is « 2.86 where R* is the radius of 
the star. The velocity scale is vq — (GM/ Rq) 1 / 2 , and the 
time-scale is to = Rq/vq. The unit of time in our plots is 
the rotational period at r — Rq: Pq — 27rfo- The other 
scaling parameters are those for the magnetic field, Bo, 
the density, po, and the pressure, pq. 

The corotation radius, r cr , is the distance at which the 




Fig. 1. — The figure shows the opening of the magnetic field lines during the disk-star interaction at T = 0, 35, 70. Time 
T is measured in Keplerian periods of rotation at r — 1 . The bottom plot represents expanded view of the magnetosphere 
near the disk at T — 70, where field lines are closed or only partially opened. 



angular velocity of the star il* equals the Keplerian an- 
gular velocity of the disk, r cr — (GAf/51, 2 ) 1 / 3 . In di- 

mensionless form, r cr = l/fi 2 / 3 , where fi* = fl^/ilo, and 
O = {GM/Rlfl 2 . We use the dimensionless angular ve- 
locity normalized to the Keplerian velocity at the surface 
of the star: w* = fi*/Q*K = (R*/r cr ) 3 / 2 . Subsequently, 



we use only dimensionless variables R — R/Rq, etc. but 
drop tildes. Results obtained in dimensionless form are 
general and can be applied to objects with widely differ- 
ent scales. 

The initial conditions are similar to RUKL. Namely, 
there is a cold dense disk with inner radius Rd and a hot 




Fig. 2. — Radial variations of the magnetic field components and the angular velocity above the equatorial plane of 
the magnetosphere of the star for = 0.58 at T = 70. In the left-hand panel the thick solid lines show B r and 
components. The thick dashed line shows the B z component. The thin solid and dashed lines show the initial dipole field 
B z and B r . In the right-hand panel the solid and dashed lines show the logarithm of the angular velocity and Keplerian 
angular velocity. 



3 






"0 2 4 8 10 12 "o 2 4^,6 8 10 12 "o 2 4^,6 8 10 12 

Fig. 3. — Episode of accretion for = 0.58. From left to right the figure shows before, during, and after the accretion 
event. The color background shows the density levels, that change from p = 0.003 (dark blue) to p = 1.6 (red). The white 
lines are magnetic field lines. The red solid line shows the surface (3 — E ma tter I 'E mag = 1 which separates the magneti- 
cally dominated region from the matter dominated region. The Alfven radius (J) = 1) is r a ~ 8. For a non-rotating star, 
r A ss 1. 



rarefied corona elsewhere. The region r < Rd rotates with 
the angular velocity of the star, while in the region r > Rd 
both the disk and the corona rotate with the angular veloc- 
ity lu = w(r) = WRr(r) so as to avoid an initial discontinuity 
of the magnetic field at the boundary between the disk and 
corona. The density and temperature distribution in the 
disk and corona are derived from the balance of the grav- 
itational, centrifugal, and pressure gradient forces. The 
disk and corona are in pressure balance. 

The initial density distribution of the disk is differ- 
ent from a stationary, non-magnetized Shakura-Sunyaev 
(1973) disk in the respect that the initial surface density 
increases gradually with distance in the outer part of the 
disk well outside of the Alfven radius. An a— type viscosity 
was incorporated in our code so as to give a controllable 
inflow of matter from the outer disk. We tried different 
a values in the range 0.01 — 0.05 and used a = 0.02 as 
the base one for most of simulations. As a result of the 
viscosity the outer disk evolves slowly towards a station- 
ary Shakura-Sunyaev disk, but there is insufficient time 
for the equilibrium to be reached. The conditions at the 
outer boundary (located at R m ax ~ 27) allow the inflow 
of additional matter to the computational region. How- 
ever, matter inflow across R m ax was not significant for the 
simulation times discussed here. 

Compared with RUKL, we gradually increased the ro- 
tation rate of the star from a small value lj* = 0.014 up 
to the final value a;* = 0.58 during several rotations of the 
star. The inner radius of the disk was placed at the corota- 
tion radius corresponding to the initial spin, Rd = r cr = 6, 
that is, far away from the star. These modifications of ini- 
tial conditions helped to overcome numerical difficulties in 
this extreme case of a rapidly rotating star. These initial 
conditions were tested and checked in multiple runs with 
different initial Rd, different times of rotational accelera- 
tion, variations of the matter distribution in the disk and 
corona, and sizes of the simulation region. 



3. RESULTS 



We performed multiple simulation runs for different an- 
gular velocities of the star a;* with the other parame- 
ters fixed. The investigated range was from cj* w 0.086 
(r cr = 1.8) (weak propeller) to = 0.82 (r cr = 0.4). 
Below we discuss a case where = 0.58 (r cr — 0.5). 

We observed, that initially the disk moved outward so 
that the inner radius of the disk increased from r = Rd = 6 
to r w 8 — 9. Then it stabilized around this value. We ob- 
served that the structure of magnetosphere around the star 
was different from the magnetospheres of slowly rotating 
stars. Namely, the magnetosphere expanded in both ra- 
dial and vertical directions. Magnetic field lines from the 
poles of the star expanded vertically and open, forming 
"magnetic towers" (Figure 1, top). Only a fraction of the 
total magnetic flux (~ 1/3) goes into the magnetic tower. 
The remainder of the flux forms a closed magnetosphere, 
which does not expand vertically, but instead expands in 
radial direction. This is shown in the lower panel of Fig- 
ure 1. The magnetic field is strongly twisted in this mag- 
netosphere so that the toroidal component is much larger 
than the other components. This is shown in the left-hand 
panel of Figure 2. The field varies as B$ ~ r -o.85 j n ^ ne 
equatorial region. Note that similar gradual distribution 
of the field was found in case of spherical accretion to the 
propelling star (Romanova et al. 2003a). 

The magnetosphere expands due to the rapid rotation 
of the star. The accreting matter stops at the Alfven ra- 
dius ta, where the matter pressure equals the magnetic 
pressure, pv^+p rs B 2 /8it . Because the B— field becomes 
non-dipolar, we can not use the standard formulae for r a , 
but instead must to find ta numerically. We found that 
the Alfven radius w 8 — 9 which is much larger than 
that for a non-rotating star [ta ~ 1)- We performed sim- 
ulations at a wide variety of and found that the Alfven 
radius increases with lu„ as rA ~ ■ 

The magnetosphere rotates with a super-Keplerian ve- 
locity, but not rigidly as shown in the right-hand panel of 
Figure 2. The angular velocity gradually decreases from 
a maximum value near the star to a Keplerian velocity of 
the disk at the magnetospheric radius ta- Thus, there is 
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no sharp transition in angular velocity between the mag- 
netosphere and the disk and "classical" propeller effect 
of pushing away incoming matter by the rigidly rotating 
magnetosphere is not observed. Instead, the disk "stands 
off" near ta and tends to be "lifted" out of the equato- 
rial plane. This leads to the formation of elongated fun- 
nel stream (see Figure 3, left panel). This stream first 
encounters the barrier formed by the rotating magnetic 
tower. However, the stream gradually forces the mag- 
netic field lines of the tower to reconnect and thus creates 
the path for accretion to the surface of the star (middle 
panel). After the accretion event, the magnetic field of 
the tower expands, forming an obstacle to further accre- 
tion (right panel). Then matter accumulates again and 
a new cycle of the accretion begins (see animation in: 
www.astro.cornell.edu/us — rus/diskjprop). This phe- 
nomenon of quasi-periodic outbursts associated with re- 
connection of the magnetic field lines was predicted by Aly 
& Kuijpers (1990) and was observed in numerical simula- 
tions of accretion to slowly rotating stars (see Goodson 
et al. 1997, 1999; RUKL; van Rekowsky & Brandenburg 
2004; Kato, Hayashi & Matsumoto 2004). 

Accretion is possible because in the magnetosphere of 
the propelling star the angular velocity decreases with the 
distance. At the distance r « ta, the centrifugal force is 
about equal to the gravitational force. This is why matter 
accretes through the funnel streams as in cases of slowly 
rotating stars. In case of propeller the variability is much 
stronger because the magnetic tower temporarily blocks 
the accretion. The quasi-period observed in the case of 
the propeller with uj* — 0.58 is about AT « 15 — 17 is 
large compared to cases of propellers with smaller angular 
velocity AT « 3 - 5 (RUKL). 

We estimated the accretion rate through the disk at the 
distance r = 10 which is beyond the Alfven radius. We 
observed that the inner regions of the disk oscillate in the 
radial direction with a quasi-period AT w 5 — 6 which 
is shorter than the main quasi-period of accretion to the 
star. The average accretion rate is M « 0.1 — 0.2. 

Angular momentum is transported between the star and 
the disk by the magnetic field of the closed magnetosphere. 



The magnetic tower influences to the dynamics of gas flow 
around the propelling star, but it is not important for the 
angular momentum transport. 

Figure 4 shows the temporal evolution of the matter 
flux to the star M, and the torques due to the matter flow 
N TO , and due to the magnetic field N/. The total torque 
is mainly due to the magnetic field. The torque due to 
the accreting matter is negligibly small even in periods of 
enhanced accretion, like in cases of slowly rotating stars 
(RUKL). We note that the star spins down continuously, 
even during periods when the accretion is blocked by the 
magnetic tower. When the accretion rate increases, the 
spin-down torque also increases. This is in opposite to the 
standard spin-down theory, where the spin-down torque 
decreases with accretion rate. The difference is due to the 
fact that during periods of enhanced accretion, a larger 
fraction of the star's flux is in the closed magnetosphere. 
This leads to a larger spin-down torque. Running a num- 
ber of cases with different o>* , we obtained the dependence 
for the average torque in the form Nf ~ — u>i . Note, that 
the same dependence was found in simulations of quasi- 
spherical accretion to the propelling star, which were done 
with different numerical code (Romanova et al. 2003a). 

The early ideas that accreting stars may spin-down due 
to the magnetic disk-star interaction were developed by 
Ghosh & Lamb (1979) and included the prediction that a 
star may accrete and spin-down at the same time. This 
last conclusion was observed in RUKL for the case of a 
weak propeller. The effect is observed in the present sim- 
ulations of much more rapidly rotating stars. The spin- 
down rate is observed to increase with accretion rate. This 
dependence is new and different from that for slowly ro- 
tating stars. 

4. DISCUSSION 

We observed that disk accretion to a star in the propeller 
regime exhibits a number of interesting new features: (1) 
The most remarkable one is the fact that the rate of spin- 
down increases as the accretion rate increases. This is 
opposite the dependence in the non-propeller regime (e.g., 
Ghosh & Lamb 1979). (2) The fact that accretion occurs in 




Fig. 4. — Temporal evolution of matter flux to the star M, and torques associated with matter N m and with magnetic 
field N f 
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the strong propeller regime is also important. The accre- 
tion is quasi-periodic and occurs through elongated fun- 
nel streams which go around the region of centrifugally- 
dominated magnctosphcrc. (3) About 1/3 of the star's 
magnetic flux goes to the magnetic tower which represents 
a region of opened field lines. However, the remainder of 
the flux goes into a radially expanded closed magneto- 
sphere, which connects the star and the disk. This mag- 
netic field is responsible for the strong angular momentum 
transport between the star and the disk. (4) The propeller 
stage may be very efficient in spinning-down rapidly rotat- 
ing accreting stars. 

It was recently suggested that the opening of the mag- 
netic field lines disconnects the star and the disk (e.g., 
Lovelace ct al. 1995) so that the spin-down of the star 
should be significantly reduced (Agapitou & Papaloizou 
2000; Matt & Pudritz 2004). Our simulations show that 
most of the magnetic flux of the star is in a closed magne- 
tosphere and is coupled to the accretion disk. The angu- 



lar momentum of the star is magnetically transported to 
the disk. Thus, the propeller mechanism may explain the 
spinning-down of the Classical T Tauri stars. 

For much stronger magnetic fields, propeller-driven out- 
flows were observed, which will be discussed in a subse- 
quent paper (Ustyugova et al. 2004). Recently, we were 
able to perform full 3D simulations of the disk accretion to 
an inclined dipole (Koldoba et al. 2002, Romanova et al. 
2003b, 2004). Next, we plan to investigate the propeller 
regime using fully three-dimensional MHD simulations. 

This work was supported in part by NASA grants 
NAG5-13220, NAG5-13060, and by NSF grant AST- 
0307817. AVK and GVU were partially supported by 
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tract MEST # 40.022.1.1.1106 and by Russian program 
"Astronomy" . Authors are grateful to anonymous referee 
for valuable comments and suggestions. 
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